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In  June,  19110,  • limited  programme  of  sky-wave  radar  measurement!  was  undertaken  to 
investigate  the  feasibility  of  measuring  saa-stste  characteristics  off  the  coast  of  Newfoundland,  The 
programme  involved  the  collection  of  data  by  the  Communications  Research  Centre's  (CRC)  Sampled 
Aperture  Receiving  Array  (SARA)  facility  at  Ottawa,  the  compilation  of  ground-truth  maps  of  wave 
height,  wind  velocity  and  ice  coverage  by  the  Centre  for  Cold-Ocean  Resources  Engineering  (C-CONE)  in 
Newfoundland,  and  data  analysis  and  interpretation  by  both  agencies.  Transmissions  were  provided,  as 
part  of  a co-operative  programme  which  involved  other  projects,  by  the  Rome  Air  Development  Centre 
transmitter  facility  at  Ava,  New  York. 

A number  of  day-time  experimental  runs  were  made  between  October,  1980,  and  April,  1982. 
Because  it  was  felt  that  ionospheric  disturbances  would  preclude  the  recording  of  useful  results,  runs 
ware  made,  as  much  as  possible,  only  during  periods  of  relatively  quiet  ionospheric  conditions. 

As  was  expected  from  the  earlier  experience  of  U.3.  workers,  maps  of  wind  direction  usually 
could  be  derived  from  the  data,  hut  maps  of  wave  heights  were  quite  sparse!  Isas  than  one-quarter  of 
the  data  ware  usaful.  In  fact,  since  the  acceptance  criterion  developed  by  the  U.8.  workers  was  found 
to  reject  virtually  all  of  the  Canadian  data;  a manual  acceptance  technique  had  to  be  employed,  jJJhila 
this  technique  has  made  possible  the  recovery  of  experimental  results,  it  would  be  impractical  tfrtvan 
operational  radar.  N. 


Determination  of  other  sea-state  characteristics,  such  as  sea-wave  spectra,  the  magnitude  of 
the  swell  component,  or  independent  measurement  of  wind  speed,  was  not  possible  from  the  data 
collected, 

Almost  all  of  the  results  wars  obtained  by  means  of  reflections  from  the  F2  layer,  with  a few 
returns  from  the  El  layer  during  a period  of  moderate  ionospheric  disturbance.  No  effort  was  made  to 
utiliae  E-reglon  reflections  since,  in  general,  these  were  not  observed  during  the  periods  of 
operation. 

A considerable  frection  of  the  enalysie  procedure!  had  to  be  directed  towarda  averaging  and 
alignment  of  the  sea-atate  spectra,  in  an  attempt  to  reduce  the  degradation  cauaed  by  ionospheric 
Doppler  offset  end  smeer.  This  problem  wee  exacerbated  by  the  necessity  for  rslatively  long 
obssrvttion  periods,  JO  to  100  seconds,  to  provlda  adsquats  resolution  in  ths  ssa-Rtate  spectra,  in 
the  presence  of  short  persistence  times  of  ionospheric  propagation  paths,  usually  only  a few  seconds 
(and  possibly  only  fractions  of  seconds). 

It  is  hoped  that,  when  the  final  enalyeia  ie  complete,  both  an  enalyaia  algorithm  and  an 
accaptance  criterion  will  emerge  which  will  indicate  the  practicability  of  operation  of  a aaa-stata 
radar  in  areas  close  to  the  auroral  cone. 


1.0  OVERVIEW  OF  THE  EXPERIMENT 


An  exploratory  programme  of  experimental  sea-stete  measurements  began  in  Canada  early  in  1980, 
involving  e number  of  agencies; 

- the  Communication#  Research  Centre  (CRC)  of  the  Department  of  Communications,  which  operated 
the  Sampled  Aperture  Receiving  Array  (SARA)  system  at  Ottawa; 

- the  Department  of  National  Dafenca,  which  provide#  funding  for  all  SARA  activities; 

- the  Centre  for  Cold  Ocean  Resources  Engineering  (C-C0RE),  at  St.  John's,  Newfoundland; 

- the  Department  of  Fisheries  and  Oceana,  which  aponaora  various  programmes  at  C-C0RE;  and 

- tha  Roma  Air  Development  Centre  (RASC)  which  provided  elgnals  from  their  transmitting 
facilities  at  Ava,  New  York. 
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In  June,  1980,  * planning  meeting  was  held  In  the  offlcaa  of  C-CORE,  attended  by 
representative!  of  each  of  the  Canadian  agenda!  and  by  Mr,  William  Randham,  from  the  University  of 
Birmingham  in  the  United  Kingdom.  Mr.  Sandham1!  addraaa  eerved  ae  an  introduction  to  the  field  of 

see-etate  reieirch  as  well  aa  a description  of  the  current  activities  at  the  University  of  Birmingham 
and  Appleton  Laboratories  [1,2, 3, 4).  The  objective  of  the  Canadian  programme,  ai  defined  at  the 
meeting,  was  to  "establish,  in  Canada,  the  technology  for  measuring  eea-state  and  other  ocean 
parameters,  using  high-frequency  skywave  radsr,  with  a view  to  improving  the  accuracy  and  timellneaa 
of  information  use  I in  aupport  of  the  operations  of  the  Department  of  National  Defence  and  the 
Department  of  Fisheries  snd  Oceans". 

A plan  of  action  was  agreed  upon,  which  encompaaseri  initial  triala  in  late  1980,  development 
of  analysia  procedures  and  conduct  of  the  experiment  throughout  1981,  and  final  analyaia  in  1982. 
Later,  thle  plan  wet  extended  to  allow  additional  observations  in  1981/82,  and  analyaia  to  continue 
into  1983,  The  enelyel*  ectivlty  is  still  in  progress, 

The  Department  of  Communications  personnel  aet  up  and  carried  out  the  experiment,  end 
converted  the  recording!  to  a data-haae  of  calibrated,  range-gated  date  stored  o..  magnetic  tape, 
C-CORE  personnel  set  up  and  operated  e treneponder  (supplied  by  CRC)  at  a site  near  St.  John1*,  and 
compiled  hlndcast  maps  of  surface-truth  data  for  each  of  the  observing  days.  Both  agendas  analysed 
and  Interpreted  the  results,  arid  will  separately  and  jointly  produce  reports  on  the  findings. 

Assistance  by  the  Rome  Air  Development  Centre,  to  provide  signal  trenemiielons , was  required 
because  the  SARA  system,  originally  built  to  carry  out  phaia-front  studies  in  support  of 
direction-finding  raiearch,  did  not  have  e high-power  transmitting  facility, 

A separata,  independent  programme  of  short-based  measurements,  using  e Coastal  Ocean  Dynamics 
Applications  Rader  (CODAR),  was  carried  out  by  C-CORE,  funded  by  the  Department  of  Fisheries  and 
Oceans. 


1.2  COVERAGE  AREA 


The  coverage  are*  for  the  experiment  it  illustrated  in  Fig.  1.  The  width  of  the  area  was 
established  by  the  beamwidth  of  the  transmitting  antanna,  approximately  65  degrees.  The  maximum  rang* 
was  limited  hy  the  signet  levels  that  could  be  achieved  and  by  the  volume  of  data  that  could  be  han- 
dled for  analytic.  In  early  trlala,  returns  were  recorded  from  the  Gulf  of  Bt,  Lawrence,  but  later 
that  area  waa  ignored  to  concentrate  effort  on  the  areas  adjacent  to  Newfoundland. 


so 


FIG.  2,  Frequency-Modulated  Continuoue-Weve 
Tranemiaelone 

The  reflected  signal,  arriving  at  the  receiver 
after  a propagation  tima  delay,  appears  to  be 
offset  In  frequency.  The  scales  ere  greatly 
exaggerated  for  clarity i the  extent  of  the 
frequency  sweep  ie  typically  50  KHsi  the  off- 
set of  the  reflected  signal  is  less  then 
60  Hs. 


FIG.  1 Coverage  Area  for  the  Canadian  Sky-Wav* 

• Sea-State  Experiment 

The  azimuthal  extent  wee  determined  by 
the  baamwidth  of  the  transmitting  antenna. 


In  the  early  trials,  problems  of  array  ambiguity  limited  non-ambiguoue  coverage  to  a narrow 
sector  near  the  northern  tip  of  Newfoundland  and  in  the  Gulf  of  Rt.  Lawranct.  A re-designed  array 
configuration  later  permitted  non-ambiguoue  recaption  over  tha  entire  illuminated  area. 


1.3  SARA  OPERATING  TECHNIQUES 

The  8ARA  system  [5]  has  tha  unique  capability  of  separately  recording  the  signal  received  by 
each  element  in  the  tntenna  array,  using  phase- locked  quadrature  detector*  in  each  rocaiver  to 
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preserve  both  phase  and  amplitude  information  for  later  conet  ruction  of  antenna  beamt  during  analy- 
•ia.  In  thie  way , the  system  can  simultaneously  obaarva  all  aaimuths  that  are  illuminated  by  the 
tranenitter. 

1.3.1  Rang*  Selection  and  Analysis 

The  ayetam  uaee  frequency-modulated  cont inuoua-wava  (FMCW)  signal*  to  permit  diacrimination  in 
range.  Tha  precisely-controlled  iwept-f requency  aignala,  after  reflection  from  a target,  return  to 
the  receiver  with  the  uaual  delay.  Aa  illustrated  in  Fig.  2,  compariton  of  tha  receive..  aignal 
againat  a locally-generated  replica  of  the  transmitted  aignal  converta  the  time  delay  to  a frequency 
offeet.  Targeta  at  different  ranges  appear  at  different  output  frequencies,  which  makes  the  extent  of 
the  obacrvable  interval  in  range  dependent  upon  the  bandwidth  of  tha  receiver. 

Diacrimination  in  range  la  accompliahed  by  tha  application  of  a complex  Fourier  tranaform  to 
tha  output  aignal,  which  dividae  tha  range  interval  into  a number  of  range  cella,  each  represented  by 
a spectral  line.  Tha  content  of  each  cell  repreaanta  the  complex  automation  of  ail  targeta,  at  that 
particular  range,  for  all  asimuths. 

Tha  bandwidth  of  tha  receiver,  and  thereby  tha  number  of  range  cella  that  can  he  observed 
simultaneously,  is  limited  by  the  maximum  sampling  rate  of  the  date-recording  system.  For  the  sea- 
state  experiment,  the  parameters  salacted  yielded  a 72  km  rango  interval,  resolvable  into  24  indepen- 
dent range  cella.  Bach  sweep  of  the  aignal  waveform  provided  a fresh  eatimata  of  the  contents  of 
thaaa  24  calls.  For  convenianea  in  the  analysis,  the  interval  was  divided  into  32  not-quite-indepen- 
dant  cells)  later  sections  of  this  paper  will  refer  to  groups  of  J2  cells, 

Tha  position  of  the  range  interval,  relative  to  tha  SARA  site,  ie  sat  at  the  time  of 
recording,  by  delaying  the  start  of  the  locally-generated  replica  of  the  sweep  by  an  appropriate 
timing  of fart,  The  signals  recorded  represent  the  contents  of  an  annulua  or  "range  ring",  about  the 
radar,  whicn  has  a radius  set  by  the  selected  value  of  range  offset  end  a range  extent  established  by 
the  receiver  bandwidth.  Of  course,  eignals  are  observed  only  from  that  sector  of  the  ring  which 
intersects  the  area  illuminated  by  the  transmitter. 

1.3.2  Beam  Forming 

Range-ring  analysis  is  carried  out  for  all  of  tha  signals  recorded  from  all  of  the  elements  in 
the  array.  From  those  results,  antenna  besais  are  then  synthesised  by  complex  sunseetion  of  the  signals 
from  corresponding  range  cella  from  each  array  clamant,  after  application  of  phase  adjustments  appro- 
priats  to  tha  selected  beam  direction.  This  procedure  is  repeated  for  each  of  the  32  range  calls, 
and,  of  course,  is  carried  out  for  all  of  tha  frequency  sweeps,  each  representing  a fresh  sample  of 
information.  The  entire  procedure  must  be  repeated  for  each  selected  beam  direction, 

The  beam  direction  is  a line  at  an  angle  relative  to  the  line  of  the  array.  When  rotated 
about  the  array  line,  the  beam  direction  describes  the  surface  of  a cone.  For  this  reason,  the  selec- 
ted beam  directions  are  commonly  referred  to  as  "cone  angles")  conversion  of  a cona  angle  to  a bearing 
on  the  earth's  surface  requires  an  estimate  of  the  elevation  angle  of  the  signal  path.  For  this 
experiment,  elevation  angles  were  estimated  externally  by  the  use  of  vertical  ionograms. 

1.3.3  Measurement  of  Doppler  Spectra 

With  reference  to  Fig.  2,  it  can  be  seen  that  Doppler  offset  of  a returning  signal  is  indis- 
tinguishable from  an  increment  in  range.  The  resulting  error  is  made  negligible  by  the  appropriate 
choice  of  operating  parameters.  However,  the  Doppler  spectrum  contains  the  information  needed  for  the 
experiment  end  so  it  mutt  be  extracted  from  the  data. 

Since  the  analyels  procedures  preserve  both  the  phaee  and  amplitude  information  of  the  signals 
representing  the  range  cells,  a complex  Fourier  transform  applied  to  a sequence  of  samples  from  a 
particular  range  cell  will  produce  a measure  of  the  rate  of  change  of  range  with  time,  which,  in  fact, 
it  the  Doppler  spentrum.  The  repetition  rate  of  the  transmitted  sweeps  and  the  duration  of  the  obser- 
vation establish  the  Doppler  bandwidth  and  the  doppler  resolution,  respectively.  Tha  transformation 
is  performed  for  all  range  cells  and  for  all  selectsd  beam  directions. 

This  final  step  in  tha  routine  analysis  generstss  a Doppler  spectrum  for  each  ocean  cell  in 
the  coverage  ares.  These  spectra  are  recorded  on  magnetic  tape  for  the  interpretive  analysis  to 
follow. 


2.0  FACILITIES 


2.1  RKCSIVIIW  BY8TIM 


The  SARA  antenna  array  comprises  90  independent  array  elements,  laid  out  in  the  shape  of  a 
cross,  as  illustrated  in  Fig.  3(a).  The  array  used  for  an  experiment  may  be  tailored  to  the 
requirements  of  that  experiment  by  the  selection  of  any  combination  of  these  elements.  Unused 
elements  usually  are  removed.  The  number  of  elements  used  in  an  experiment  is  limited  ultimately  by 
the  maximum  aampltng  rate  of  the  data-racording  system,  taking  into  consideration  a number  of 
interrelated  factors  such  as  resolution  in  range  and  Doppler,  rang*  window,  Doppler  bandwidth, 
acceptable  levels  of  aliasing  in  range,  and  tha  available  selection  of  sweep  rates  and  receiver  output 
filter*.  For  the  sea-state  measurements,  the  parameters  chosen  dictated  an  array  sis*  of  42  elements. 

Two  42-*l*m*nt  configurations  were  tried,  At  first,  in  an  attempt  to  obtain  the  minimum 
possible  beamwidth  (less  than  1 degree),  a spars*  selection  of  elements  was  made  along  the  full  extent 
of  the  array,  a*  illustrated  in  Fig,  3(b).  Tha  32  closely-spaced  central  elasMnt*  were  combined  in 
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pain,  reducing  Chair  numbar  Co  16  to  allow  all  26  of  tha  outer  alananta  Co  ba  uaad,  Although  thia 
configuration  waa  known  to  ba  ambiguoua,  it  waa  thought  that  tha  amblgultlea  could  ba  raaolved  or 
ignorad,  at  had  haan  poatible  in  direction-finding  axparimanta.  To  a large  extent,  thia  waa  not  tha 
cate,  however,  and  ao  tha  array  had  to  ba  reconfigured.  Soma  reaulta  ware  reacued  from  thia  early 
attempt  by  deleting  the  aignala  from  the  outer  elementa  during  the  analyait  procedurat,  producing  a 
wider  beam  pattern  with  a capability  for  non-ambiguoua  reception  from  a email  aector  in  the  middle  of 
the  intended  coverage  area. 


L'V 

i 


Although  aignala  ware  received  from  tha  entire  array,  ambiguity  problama  required 
that  only  the  central  portion  of  the  array  be  uaed  for  analyaia. 

I 

i 


v 


Tha  elaaanta  of  tha  cloaaly-apaced  centra  portion  were  uaad  Individually  for  tha 
final  configuration. 
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The  dealgn  of  the  Mcond  configuration,  illustrated  in  Fig.  1(e),  took  that*  probltai  into 
account  and  accepted  the  wider  baamvldth  of'. a ahortar  array  to  avoid  aabiguoul  raaponaaa  within  tha 
covaraga  araa.  To  raduca  aidaloba  raaponaaa,  tha  deeigr  of  tha  array  naeaaaarily  included  tha 
application  of  an  antanna-ihading  "window"  function.  In  thia  caaa  a coaina-aquarad  window,  appliad  aa 
a function  of  diatanca  along  tha  array,  auppraaaad  aidaloba  raaponaaa  to  a laval  about  20  decibel* 
balow  tha  raaponaa  of  tha  main  loba.  Tha  baamwidth  of  tha  ahortanad  array,  with  tha  "window"  applied, 
waa  about  2 dagreaa. 

For  analyala  of  tha  limited  covaraga  aactor  of  tha  original  configuration,  tha  tana  window 
function  waa  applied  to  the  aame  array  extant  that  waa  uaad  in  tha  aacond  configuration,  and  tha 
raaponaaa  of  tha  outlying  aleioenta  ware  aat  to  aaro.  Thia  produced  a haam  pattern  nearly  identical  to 
that  of  the  aacond  configuration,  except  for  tha  anbiguoua  raaponaaa  that  emerged  whan  the  beam  waa 
ataarad  outalda  tha  limitad  covaraga  araa. 

Tha  SARA  racaiving  and  data-racording  ayatam  ii  illustrated  in  Fig,  A,  Important  faaturea  of 
tha  ayatam  ara: 

(1)  tha  antenna  aalaction  panel,  where  tha  aalactad  array  alementa  ara  interconnected  with  tha 
bank  of  racaiverai 

(ii)  tha  computer-controlled  local-oacillator  and  fraquancy-awaep  generator,  which  perform  tha 
function*  of  timing  tha  racaivara  and  da-ramping  tha  awtpt-fraquancy  aignala  by  providing 
phaaa-matchad  local-oacillator  and  rafaranca  aignala  to  all  racaivara; 

(iii)  tha  aampla-anri-hold  unit  which  fraaiaa  tha  valuaa  of  all  of  tha  receiver  output  aignala  whan 
triggarad  by  a aampling  pulaa  from  the  timer; 

(iv)  tha  multiplexer,  which  aaquantially  aeana  the  voltagaa  held  by  tha  eample-and-hold  unit, 
following  each  aampling  pulaa; 

(v)  tha  analogua-to-digital  converter,  which  convarta  tha  aignala  to  12-bit  digital 
rapreaentation; 

(vi)  tha  computer,  which  managea  tha  taek  by  controlling  timing,  avant  aaquenca  and  frequency 
aalaction,  and  formata  the  data  into  racorda  for  transmission  to  tha  tape  recorder; 

(vll)  tha  frequency  atandard,  to  which  all  frequency-generation  an^  timing  functiona  ara 
aynchronlxed. 


FIG.  A.  Tha  SABA  he ca lying  and  Date-Recording  8yatam 

Tha  required  high  pracialon  of  ilgnal  frequency  aalaction  and  timing  control  la 
aatabliahcd  by  the  rubidium  frequency  standard. 


Tha  racaivara  ara  locaUy-deoigned  double-iuparheterodyne  unite,  tuned  aeroaa  tha  2-30  HHa 
band  by  computer  aalaction  of  local-oacillator  fraquancisa.  Each  receiver  containa  a quadrature 
detector  which  parforma  a third  frequency  translation  to  ahift  the  output  centra  frequency  to  laro, 
with  "poaltive"  and  "nagativa"  frequency  band*  resolvable  from  tha  quadrature  aignala.  Tha  output 
bandwidth  of  aach  racaivar,  aatablliliad  by  intarchangaabla  poat-dataotion  low-paai  filter*,  i*  quite 
narrow.  For  tha  aaa-atata  axparlmant,  tha  bandwidth  uaad  wa*  +60-H*.  Sach  racaivar  haa  an  indepen- 
dent automatic-gain-control  (AGC)  circuit.  To  enable  restoration  of  tha  input  signal  voltagaa  during 
analysis,  tha  individual  AGC  control  voltagaa  were  alao  campled  and  recorded. 
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Th#  aignala  from  aach  receiver  had  to  ba  aaaplad  at  a rata  high  enough  to  avoid  algnificant 
aliaaing  ("Nyqulet  folding").  For  tha  aaa-atata  experiment  tha  two  output  aignala  and  tha  AQC  voltaga 
from  aach  of  tha  42  vacaivara  had  to  ba  sampled  at  a rata  conmanaurata  with  tha  60-Ha  filtara.  Thla 
worked  out  to  a total  sampling  rata  of  25200  aamplts  par  aeeond.  Thia  waa  aaally  handlad  by  the 
analogua-to-digital  convartar,  which  can  oparata  at  rataa  up  to  45000  aamplaa  par  eaconri,  but 
approached  tha  Uniting  rata  of  tha  tapa  racordar,  which  can  accapt  up  to  about  27000  data  worda  par 
aacond. 

lonograma  ware  received  by  the  aama  ayaten,  using  a aaparata  receiver  and  a aaparata  antenna. 
The  output  bandwidth  of  the  ionogram  receiver  ia  wide,  to  permit  obaervation  of  the  entire  range  of 
interest,  lonograma  were  analysed  on-line,  for  immediate  presentation,  and  were  aiao  recorded  to 
permit  later  re-analyala  if  neceasary. 

Tha  Honeywell  DDP516  computer  ayatan  was  alio  uaad  for  analyaia  taaka,  both  on-line  and 
off-line.  It  ia  quite  fait,  daapita  ita  age,  and  ia  limited  in  capability  only  by  tha  relatively 
small  capacities  of  ita  memory  and  diac  racordar,  and  by  tha  availability  of  only  ona  tapa  racordar. 
Tha  computing  facilities  include  a dot-addressable  dot-matrix  printer  which,  in  concert  with  a full 
ranga  of  graphics  software,  can  ganarata  a variety  of  graphic  prlnt-outa,  including  Doppler  spectra 
and  gray-ecale  lonograma,  all  at  high  spaada. 

2.2  TRANSMITTER 

Tha  target  area  in  tha  octan  waa  illuminattd  by  aignala  tranamittad  from  tha  RADC  tranamlttar 
facility  at  Ava,  Haw  York.  That  facility  haa  available  both  a aalaction  of  transmitters  and  a 
aalactlon  of  antannaa.  RADC  ia  llcansad  for  tha  transmiasion  of  low-powar  continuous  awaapa  across 
tha  ant  Ira  band,  and  for  high-powar  narrow-band  awaapa  within  a aariaa  of  200-KHa  bands  diatributad  at 
approximately  2-HHa  intervals  across  tha  band.  Tha  low-powar  continuous  awaap  mods  ia  used  for  tha 
generation  of  lonograma,  tha  high-powar  mods  for  radar  opsrationa, 

The  aaa-atata  axparimant  uaad  a rotatable  horlsontally-polarieed  log-periodic  antanna  which 
had  an  aiimuthal  bearawidth  of  about  &3‘  ovar  ita  operating  range  of  6.5  to  30  MHi.  Thia  antanna  waa 
mounted  on  a 35-mstre  tower,  which  raaultad  in  aoma  unavoidable  areas  of  poor  illumination  btcauie  of 
nulls  in  tha  elevation  pattern.  Power  lavala  ware  limited  to  3-10  KW  for  lonograma j 10-13  KW  for 
radar.  Although  mors  radar  power  would  have  improved  tha  aignal-to-noiaa  ratio  in  aoma  casta,  thia 
waa  not  possible  bacauaa  high  power  caused  flash-over  and  damage  to  the  antanna  at  tha  higher 
frequencies. 

2.3  TRANSPONDER 

A tranaponder,  built  by  Stanford  Raaaarch  International,  waa  installed  at  8t.  John's, 
Newfoundland,  to  provide  a calibration  point  in  hearing  and  ranga,  Tha  unit  waa  oparatad  continuously 
on  all  observing  days,  retransmitting  all  racaivad  aignala  with  8.3-Ha,  supprassed-carrtar  amplitude 
modulation  to  ganarata  a falsa-ranga,  falsa- Doppler  signature  in  tha  Doppler  apactra  whenever  tha 
analyaia  procedure*  selected  tha  range  call  aurrounding  St.  John'a.  A directional  array  of  Reverage 
•lament*  waa  uaad  ai  tha  antanna  to  provide  a largt  crosa-atctional  area  in  tha  direction  of  Ottawa. 

3.0  ANALY8I8  TECHNIQUES 

Analyaia  of  the  data  recorded  in  tha  axparimant  waa  aeparatad  into  four  taaka i 

U)  Analyaia  and  presentation  of  lonograma,  from  which  operating  pointa  ware  aalectad  and  virtual 
heights  wars  derived; 

(ii)  Convaraion  of  tha  recorded  narrow-band  "chirp"  radar  aignala  to  a data-base  of  calibrated 
ranga  ringa; 

(ill)  Baam-formiiig  and  dopplar  analyaia  of  tha  range-ring  data  to  dafina  Doppler  apactra  from  indivi 
individual  ranga  calla;  and 

(iv)  Interpretation  of  Dopplar  apactra  to  derive  sea  conditions. 

The  first  two  of  thsie  tasks  wars  carried  out  antirtly  at  CRC/SARA;  the  latter  two  wara 
undertaken  independently  by  both  CRC  and  C-CORE.  CRC  effort  waa  concentrated  upon  technique*  of 
maaiuring  wave-height;  C-CORE  effort  upon  technique!  of  automatic  generation  of  wind-field  maps  by 
maasuramant  of  wav*  direction. 


3 . 1 Analyaia  and  Presentation  of  Ionogr«ms 

lonograma  war*  analysed  and  praaanttd  in  real  time  by  the  flARA  control  computer,  Although  not 
equipped  with  a hardware  fast-Fouriar  procasaor,  the  Honeywell  computer  ia  fait  enough  to  convert 
Incoming  awapt-fraquancy  aignala  to  lonograma  in  real  time.  Fig.  5 i*  a reproduction  of  a gray-acal* 
backacattar  Ionogram  which  waa  drawn  originally  on  a plain  paper  ahaat  approximately  28  X 30  cm,  uaing 
th*  dot-addraasabi*  matrix  printer  attached  to  the  computer. 

A new  ionogram  waa  generated  every  13  minutes,  followed  aach  time  by  tha  recording  of  two 
radar  runs  at  two  ranges  aalectad  from  th*  ensemble  that  covered  th*  area  of  interact.  Both  th* 
ionogram  aignala  and  th*  radar  aignala  war*  ratalnad  on  tapa. 
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FIG , 5 1 Backecatter  Ionogram 

Ihaaa  ware  generated  avary  fiftaan  minute*  during  experimental  run*.  The  SARA 
control  computer  performed  the  required  Fourier  analyala  and  converaion  to  a 
graphic  print-out  in  real  time, 


A portion  of  an  overhead,  nearly-vart leal  ionogram  appeared  in  tha  lower  laft  corner  of  each 
backecatter  ionogram,  raaulting  from  illumination  of  tha  overhead  ionoiphar*  by  aidalob*  radiation 
from  tha  tranamitting  antenna,  Only  the  upper  and  of  tha  overhead  trace  waa  available  bacauae  of  tha 
6.5  MHa  lower  operating  limit  of  tha  tranemitting  antenna,  but  thia  waa  uaually  enough  to  determine 
the  virtual  height  of  tha  overhead  lonoaphare  at  the  aalected  operating  point,  taking  into 
eonaideratlon  tha  uaual  tranaformation  for  angle  of  incidence.  Aa  1*  often  the  ceae  in  aaa-atate 
experiment*,  the  ionoaphera  above  the  receiver  waa  aaaumed  to  be  a valid  indicator  of  ionoapherlc 
condition*  at  the  actual  radar  reflection  point  about  1000  Km  away, 

The  backacattar  trace,  extending  outward*  from  the  doubla-hop  overhead  trace,  wee  ueed  for 
aelectlon  of  the  operating  point.  A technique  for  the  optimum  aalection  of  algnal  frequency,  baaed  on 
the  uee  of  a nomogram  overlay  in  conjunction  with  a complete  overhead  ionogram,  hae  bean  deecribed  by 

U.8.  worker*  [6],  Although  tha  required  vertical  ionogrem  could  have  baan  generated  by  tha  uaa  of 

CRC/8ARA  tranemitting  equipment,  uaa  of  that  technique  would  have  provided  little  actual  benefit, 
becauee  of  the  vary  limited  llet  of  available  operating  frequenciea.  For  thia  reaeon,  operating 
point*  were  choaen  atmply  by  aalecting  the  higheat  available  operating  frequency  that  placed  the  rang* 
of  lntereat  behind  the  leading  edge  of  the  backacatter  trace. 

3.2  Generation  of  the  Range-Ring  Data  8aae 

Converaion  of  the  recorded  algnala  to  a ant  of  calibrated  rang*  ring*  waa  the  routine, 
although  tima-coneumlng,  taek  of  Fourier  tranaformation  of  all  of  the  de-rampad  narrow-band  aweepa 

from  all  of  the  receiver*.  Thia  procea*  converted  the  312  individual  "chirp"  aweepa,  recorded  during 

each  102.4  aecond  obaervation  dwell,  to  312  time  (ample*  for  each  of  tha  32  rang*  celle.  The 
raaulting  range-ring  data  file  for  each  radar  obaervation  comprised  312  complex  (ample*  of  32  range 
cell*  for  each  of  the  42  antenna  element*. 

The  recording  rata  waa  high  enough  to  permit  a eubetantial  "guard  band"  ( 7 1 for  protection 
againat  Nyquiet  folding)  algnala  aliaaad  from  outeid*  the  accepted  rang*  interval  were  euppreaeed  more 
than  23  decibel*.  The  data  baea  waa  compiled  at  the  CRC  computing  centre  at  a rat*  of  about  A to  S 
data  tape*  per  week,  requiring  a period  of  aevaral  week*  to  proceet  the  recording*  froea  on*  obaervlng 
day. 

3.3  BKAM-FOkHMO  AND  DOPPLER  AHALTBIB ■ PRO 

At  CRC,  the  taek  of  beam-forming  and  Doppler  analyala  waa  performed  aeperataly  from 
interpretation  in  order  that  proceduroa  requiring  large  atorag*  array*  could  be  carried  out  at  the  CRC 
computing  centre,  while  procedure*  involving  apectral  (election  and  alignment  could  be  handled  on-lina 
at  the  8ARA  computer  conaol*. 
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3.3.1  Be ta-Forming 

Antenna  beam!  were  formed  by  th«  complex  eummetlon  of  algnala  sampled  at  the  various  elements 
along  the  array,  after  appropriate  adjuitmants  in  phaae.  Use  of  this  technique,  rather  than  the 
equivalent  but  more  efficient  technique  of  Fourier  tram  format  ion,  wae  required  becauee  the  elements 
in  the  array  were  not  uniformly  spaced.  An  advantage  of  the  complex-arithmetic  technique  waa  the 
capability  for  arbitrary  aolection  of  eslmuthal  directions,  independent  of  the  array  direction  or  the 
radio  frequency. 

Beams  were  formed  for  a iixed  array  of  directions,  at  two-degree  intervale,  from  a bearing  of 
51  degrees,  which  is  part  way  up  the  Labrador  coeat,  to  81  degrees,  which  is  south  of  Newfoundland. 

The  southern  limit  was  sometimes  extended  to  a bearing  of  93  degrees,  which  wet  the  limit  of  the 
illuminated  area, 

3.3.2  Doppler  Analysis 

Using  fast-Fouriar  techniques,  Doppler  spectra  were  derived  for  each  range  cell  in  tech 
selected  beam  diractlon,  and  for  each  obsarvation  intarval. 

Although  the  length  of  the  obaervetion  intarval  had  been  aet.  at  102.4  aeoonda  at  ricording 
time,  it  wae  known  that  there  would  be  problems  of  interpretation  of  a non-a tat  ionary  procaaa  bacauae 
individual  ionoepharic  paths  frequently  exlat  for  only  a few  eeconds  at  a time,  particularly  at  high 
magnetic  latltudae  [8].  Fourier  enalyela  implicitly  assumes  a statistically  stationary  eigntl,  but  a 
long  obsarvation  of  a aky-wave  signal  may  encompass  the  sum  total  of  several  propagation  conditions 
which  have  appeared  end  disappeared  at  different  times  during  the  Interval,  in  the  finel  result, 
these  appear  ee  multipath  components.  To  snelyss  the  rucorded  dsta,  a compromise  had  to  bs  struck 
betwsan  t relatively  long  averaging  time,  which  provided  fine  resolution  In  the  Doppler  apactra,  and 
dagradar.ion  of  ths  spectra  raaulting  from  the  non-ntatlonary  process.  All  of  the  CRC  analysis  has 
been  carried  out  on  half-intervals  of  tha  recorded  data,  l.a.,  each  102.4-aacond  recorded  interval  has 
been  treated  at  two  contiguous  51.2-second  intervale. 

To  eupprett  spectral  leakage,  with  a minimum  of  lint  broadening,  e Rltckman-Herrie  minimum 
three-term  window  [9]  was  applied  to  each  of  tha  half-interval  segment*  of  the  data.  Thie  brought 
about  the  required  improvement  but  reeulted  in  the  auppraaaion  and  effective  loaa  of  a considerable 
amount  of  data  at  the  ends  and  In  tha  middle  of  the  observation  interval.  By  re-analysle  of  the 
half-interval  between  the  one-quarter  and  thraa-quarters  points  of  the  interval,  and  application  of 
the  tame  window,  a third  time  Intarval  waa  defined,  permitting  the  recovery  of  much  of  tho  lost  data 
and  the  generation  of  a third  spectrum. 

The  three  apactra,  representing  time  intervale  which  overlapped  SOX,  actually  were  only  10X 
correlated  because  of  the  effect  of  the  window  function.  Hence  each  rang*  cell  wee  represented  by 
three  nearly-lndapandant  Doppler  spectra,  producing  a total  of  96  epectra  for  the  32  range  cells  in 
each  eslmuthal  direction,  for  eech  observation  interval.  These  were  recorded  on  e tap*  readable  by 
tha  SARA  control  computer. 

3.4  INTERPRETATION  OF  DOPPLER  8PECTRA,  CRC 

The  Doppler  ipectrum  of  a rider  signal  reflected  from  ocean  waves  i*  characterised  by  two 
"Bregg"  peaks:  strong  spectral  llnaa  raaulting  from  focuaeing  of  the  reflected  energy  by  the  regular 

pattarn  of  the  wavta,  in  a mechanism  analogous  to  a diffraction  grating.  Theae  peak*  are  superimposed 
upon  a pedestal  of  aacond-ordar  reflections  which  form  a continuum  of  raturna  across  a spectral  width 
of  about  2 Ht.  In  theory,  several  measure*  of  sea  conditions  can  be  derived  from  this  Doppler 
signature,  Theae  include  wav*  height,  wave  direction,  wind  spttd,  tea-wave  apectrum,  dominant  wave 
period,  and  velocity  of  ocean  currents  [10].  Since  alt  of  these  are  related  to  soma  aspect  of  the 
Doppler  signature,  the  accuracy  of  measurements  made  by  a aky-wave  radar  la  degraded  by  the  varitblt 
Doppler  offset  and  Doppler  ameer  imposed  upon  the  signal  returns  by  motion  of  the  Ionosphere.  In 
particular,  from  signals  raflecttd  by  tha  ionoepharic  F layer,  only  two  tea  conditions  may  be 
observable  with  a useful  degree  of  confidence:  wave  height,  derived  from  the  ratio  of  the  power  of 

the  prominent  Hragg  lint  to  the  power  of  tha  surrounding  continuum;  and  wav*  (wind)  direction,  derived 
from  the  ratio  of  the  power  of  the  "approaching"  Bragg  line  to  that  of  the  "receding"  Bragg  line. 

A greet  deal  of  subjective  analysis  it  necessary.  To  btgin  with,  tha  redaction  process  at 
the  sea  surface  displays  considerable  variance,  necessitating  tha  incoherent  averaging  of  many 
individual  apactra  derived  from  several  independent,  cloaely-epaccd  time  sample*  from  several 
Independent,  cloaely-ipaced  range  ctlle.  In  the  absence  of  ionospheric  effect*,  simple  averaging  of 
group*  of  up  to  100  spectra  usually  can  resolve  all  of  tha  at*  conditions  listed  above,  However,  tho 
ionospheric  Doppler  anomalies  imposed  upon  aky-wave  radar  raturna,  if  unconditionally  averaged  over  a 
large  number  of  apactra,  tend  to  degrade  the  final  result  to  a smeared  signature  from  which  little  or 
nothing  can  be  derived.  Ionospheric  Doppler  offeet  and  the  incidence  of  multipath  propagation  both 
vary  considerably  ai  functions  of  ttme  and  space,  Their  effects  can  be  minimised  by  adjustment!  of 
the  unconditionally-averaged  spectra  prior  to  final  averaging:  spectra  displaying  significant 

multipath  tffacte  can  be  discarded;  those  showing  Doppler  offsets  can  ba  shifted  to  align  the  Bragg 
line*  of  the  sea-tcho  epectra.  Several  observation:!  of  tha  seme  call  are  frequently  necessary  to 
obtain  a useful  result, 

Tha  high  varlenc*  of  the  reflection  and  propagation  process**  imposed  a dilemma:  the  averaged 

results  were  needed  to  accurataly  assess  the  ionospheric  affects,  but  those  effect*  had  to  be 
minimised  to  obtain  the  average,  Tha  working  compromise  employed  unconditional  averaging  of  subsets 
of  the  ensemble  of  individual  spectra,  to  produce  a aet  of  epectra  useful  for  estimating  ionospheric 


effects,  followed  by  selection,  alignment  and  conditioning  of  th«  unconditional  avaraRaa  prior  to  tha 
calculation  of  the  final  average  fron  which  aea-etata  charactarlatlca  were  derived. 

At  CRC  tha  96  apactra  from  an  obaarvation  Interval  ware  flrat  unconditionally  averaged  in 
groupa  of  12:  three  "overlapped"  apactra  from  each  of  four  contiguoua  range  cells,  Thia  produced 

eight  unconditionally-averaged  apactra  fron  each  obaarvation  interval.  Theae  were  examined 
individually  by  the  experimenter  on  the  conaole  oscilloscope  of  tha  SARA  computer,  aaelited  by  a 
controllable  flashing  display  which  indicated  the  expected  separation  of  tha  Bragg  lines.  Each  of  the 
eight  spectra  waa  either  rejected  or  tagged  to  indicate  the  experimenter's  best  estimate  of  tha 
location  of  tha  dominant  Bragg  line.  The  final  average  spectrum,  assumed  to  represent  conditions  for 
the  entire  32-celt  range  Interval,  waa  derived  from  tha  accepted  spectra  after  alignment  of  the  Bragg 
lines  according  to  the  experimenter's  tags,  and  normalisation  in  amplitude  in  accordance  with  the 
geometric  mean  of  the  power  in  the  dominant  Bragg  lines.  The  process  is  illustrated  in  Riga.  6(a), 
6(b)  and  6(c). 

Soma  poat-averaging  selection  waa  also  done;  in  particular  if  fewer  than  three  unconditional 
averages  were  Included  in  tha  final  average  the  measurement  of  wave  height  waa  not  considered  valid. 

Before  this  manual  taehniqua  waa  adoptad,  conaidarable  effort  was  expended  upon  attempt!  to 
machanlaa  tha  aalaotion/al ignmant  process  but  this  met  with  very  little  aucceaa.  The  definition  of  en 
acceptable  epectrun  wae  difficult  beceuee  of  the  wide  range  of  acceptable  epectral  ehapee,  end 
automatic  method*  of  determining  the  poaition  of  the  dominant  Bragg  line  frequently  (elected  incorrect 
peaks.  The  technique  of  rejecting  apactra  on  tha  baala  of  the  "aquivalant  width"  of  the  dominant 
Bragg  Una,  used  by  tbe  U.S.  workara  (11)  in  their  automated  echame,  tended  to  reject  all  of  tha 
apactra  recorded  by  this  experiment,  Although  the  manual  tachnique  is  cumbersome  and  tim*-eoniumin|, 
it  hag  managed  to  extract  useful  results  from  relatively  poor  data  and  la  providing  both  statistics 
and  experience  which  may  laad  to  tha  design  of  an  automatic  procadurt. 

Wav*  heights  reported  in  thie  paper  ere  thoie  defined  ae  "elgnificent  wav*  height*",  l.e. 
r.m.e.  wav*  height!  multiplied  by  four.  They  were  eetiiaatad  by  use  of  the  Mareeca-Oeorgee  power-law 
relation  t 12] 


*Rb 

h ■ , hk0  >0.2 

hn 


where  h ■ r.m.e.  wave  height, 

R * unweighted  ratio  uf  second-order  to  first-order  power  encompassing  the  dominant  Bragg 
line, 

k0  ■ radar  wavenumber, 

a,b  ■ constant!,  0.A  and  0.6,  respectively,  found  by  Hareeca  and  Oaorgae  to  ba  tha  beat 
overall  power-law  fit  to  their  data, 

Use  of  this  taehniqua  implies  a knowledge  of  the  location  of  tha  nulls  which  separate  the 
Bragg  line  from  tha  second-order  continuum,  but  line*  theaa  null*  are  raraly  observed  In  aky-wavn  date 
(naver  in  tha  Canadian  data)  they  were  assumed  to  lie  at  0.07  Ha  (13)  each  side  of  tha  peak  of  the 
Bragg  line. 

Wavt  directions  were  aetimeted  from  the  ratio  of  the  power  in  the  two  Bragg  line*  by  uae  of 
the  Iong-Trlsn*  [14]  formula! 

(0.56  ♦ 0,5  cos  28  \ 

— - — ) ♦ 34.02  decibels 

’ ) 


where  p » Bragg  Una  power  ratio 

6 - angle  between  the  radio-wave  propagation  direction  and  tha  mean  eaa-wav*  direction. 


3 , 5 AUTOMATIC  AHALYBI8  TECHNIQUES , C-CORE 

At  C-C0RB,  analyali  effort  we*  directed  towards  the  development  of  automatic  methods  for  tha 
extraction  of  wind  (wave)  direction.  Tha  analysis  system  derive*  wind  direction*  for  a specific  eat 
of  geographical  locations,  for  comparison  against  hindcaat  maps  of  wind  directions  obtained  frost 
meteorological  data. 

initially,  attempt*  war*  mada  to  devise  tn  automatic  schema  to  maasur*  both  wav*  height  and 
wav*  direction  but,  for  tha  data  avail able,  it  was  found  that  tha  algorithm*  could  not  consistently 
identify  the  epectral  chtracterietlce  required  for  the  meaiurement  of  wav*  height,  Nearly  ell  of  the 
wav*  height  reeulte  had  to  be  rejected,  However,  for  the  meaeurament  of  wave  direction,  eufflcient 
success  was  achieved  to  warrant  continued  anelyels. 


Thu  final  spectrum  resulting  from  the  processing  of 
the  four  aelacted  ipectre  from  Fig.  6(b)  he*  been 
centred  for  presentation.  Analysis  of  Its  characterls- 
tice  yielded  a wave  direction  of  23.4  dagreaa  end  e 
wave  height  of  3.3  metres!  both  of  which  agreed 
closely  with  hindcest  surface  data. 
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Th*  analysis  procedure  divided  Che  102.4-aacond  obtervetion  intervel,  repreeenced  by  S 1 2 
complex  samples  from  eech  receiver,  into  seven  sub-intervals  of  128  complex  samples  each.  The  sub- 
intervals  were  502  overlapped,  Eech  eub-lntervel  vae  transformed  to  a spectrum  after  application  of  a 
BlaclMan-Harria  window  function.  The  eevsn  apectra  from  each  range  cell  were  unconditionally 
averagedj  then  the  averaged  apectra  from  each  of  eight  contiguous  range  cells  wars  aligned  and 
averaged  to  create  one  final  spectrum  representing  one-quarter  of  the  72  Km  range  interval.  Using  the 
Long-Tritna  formula,  a wind  direction  was  calculated  from  the  ratio  of  the  two  Bragg  peaks,  The 
apectra,  and  the  resulting  wind  maps,  were  automatically  plotted  for  examination  and  post-analysis 
acceptance . 

Attempts  made  to  Introduce  a quality-control  index  have  been  unsuccesaful.  An  index  bated 
upon  the  normal ited  ratio  of  adjacent  peaks  and  troughs  within  the  spectrum  has  been  developed  as  a 
possible  indicator  of  ionoapherie  contamination,  but  it  remains  to  be  teen  whether  it  can  be  used  aa  a 
satiafactory  poat-analyaia  acceptance  criterion.  It  is  more  tolerant  than  the  Oeorges-Marssca 
"equivalent  width"  111)  technique,  which  rejects  virtually  all  of  the  Canadian  data,  because  of  the 
lack  of  an  acceptable  mechanism  for  rejection  of  poor  results,  quality  control  hat  bean  carriad  out  by 
visual  inspection  of  the  plotted  apectra,  after  analyaia.  Hind  directions  ware  retained  only  for 
accepted  spectra. 

Automatic  methoda  of  alignment  of  the  unconditionally  averaged  apectra,  aa  a meant  of  reducing 
ionospheric  effects,  still  pose  problems,  Attempts  were  made  to  align  apectra  by  determination  of  the 
position  of  the  dominant  Bragg  peak,  taking  into  account  the  content  of  the  target  cell  (land,  coast- 
line or  ocean),  but  groaa  misalignments  tended  to  occur  if  noisy  apectra  wars  intermixed  with  good 
spectra  or  it  adjacent  apectra  represented  different  types  of  targats. 

The  alignment  technique  employed  made  use  of  measures  of  the  average  value  and  the  centroid  of 
each  spectrum  to  estimate  the  required  shift  in  frequency.  Beginning  at  the  centroid,  the  algorithm 
searched  outwards,  in  each  direction,  until  it  found  a point  where  the  power  dropped  below  the  average 
value.  The  midpoint  between  tha  limits  thus  found  was  assumed  to  be  the  correct  midpoint  of  tha  spec- 
trum for  purposes  of  alignment.  The  advantage  of  tha  use  of  this  method  was  its  relative  insensiti- 
vity to  the  effects  of  target  types  and  mixtures  of  target  types,  e.g.,  along  a coastline  one  call 
could  be  on  land  and  the  next  on  the  ocean.  The  technique  caused  a certain  amount  of  smearing  of  the 
Bragg  linea,  but  appeared  to  be  suited  to  the  alignment  of  the  short  (128  point)  intervals  used  in  the 
analyses, 

CRC  results  indicate  that  satisfactory  wind  directions  frequently  cart  be  derived  from 
unconditional  apactral  averages;  this  suggests  that  a simplification  of  the  C-CORE  techniques  might  be 
possible. 

4.0  RESULTS , CRC 

CRC  analysis  effort  was  focussed  upon  the  derivation  of  wave  height.  Of  course,  any  spectrum 
that  yielded  an  acceptable  measure  of  wava  height  also  yielded  a wave  direction.  In  fact,  it  was 
possible  to  make  meeaureatenta  of  wave  direction  even  in  the  presence  of  severe  multipath  distortion 
because  the  analysis  technique  permitted  the  experimenter  to  atipulate  which  spectral  lines  were  the 
Bragg  peeks.  Aa  a result,  the  CRC  analysis  procedures  produced  complete  wave-direction  maps  ae  i 
by-product  of  tha  wave-height  enalyaii.  On  the  othar  hand,  tha  wave-height  mips  displayed  only  a 
sparse  distribution  of  results  because  of  the  difficulty  of  obtaining  sufficient  numbers  of  acceptable 
unconditionally-averaged  spectra  to  generate  e useful  final  spectrum. 

Two  observing  daye  were  analysed:  April  9,  1981,  when  waves  vara  generally  about  two  maLrea 

in  height;  and  April  8,  1982,  whan  an  Atlantic  atom  south  of  Nowfoundland  produced  wave  halghta  of  4 
to  7 metres  within  tha  radar  coverage  area. 

It  should  be  noted  that  tha  aucetas  to  be  daicribed  In  the  enelyaie  of  theee  two  selected  run* 
was  not  generally  attainable.  Preliminary  analyaia  of  othar  runs  indicated  that  many  of  them  will 
produce  little  or  no  results. 

4.1  April  9.  1981 

Hava  height*  measured  on  April  9,  1981,  are  shown  in  Fig.  7,  overlaid  on  conventional  hlndcaat 
wava  height  contours  derived  from  meteorological  data  and  ship  reports.  The  radar  results  wars 
obtained  using  tha  ambiguous  array  configuration,  and  to  the  southern  half  of  th*  plotted  data  aey 
have  been  Influenced  by  signal*  received  from  an  area  near  th*  northern  edge  of  th*  map.  Tha  northern 
half  of  the  plottod  data  la  free  from  ambiguous  returns. 

Tha  mtaaurad  wava  heights  wara  somewhat  high,  varying  from  correct  values  to  about  double  tha 
correct  values.  Bta*  toward  high  readings  it  to  be  axpaetad  if  powar  from  the  focussed  Bragg  lines  it 
shifted  into  tha  aecond-ordtr  continuum  by  tha  affects  of  ionospheric  multipath  propagation.  This  it 
particularly  Important  whan  wave  halghta  are  low,  because  tha  signal  power  returned  in  th* 
second-order  continuum  l*  then  relatively  low,  making  tha  ratio  of  Bragg  lint  power  to  continuum  powar 
very  sensitive  to  contamination. 


I 


Tha  radar  aaaaunmanta  art  both  nalf-conalatant  and  conalatent  with  tha  netaorolojleal 
laobara. 
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The  corraapondlng  wave-direction  map  ti  ehown  in  Fig,  8,  overlaid  on  hindcaet  preaaure 
ieobere.  Sal f-coneieCancy  of  the  reaulta  ia  clear,  ae  ia  conaiatency  with  the  air  preaaure  contoura. 
Wave  direction  mapa  produced  by  a alngle  radar  neceaaarily  diaplay  a right/left  ambiguity,  but  alnce 
thia  waa.  rcaolved  by  the  weather  map  it  waa  not  ehown  here. 

Three  -baervation  paaaea  were  made  over  the  coverage  area.  Only  the  reaulta  from  the  aecond 
paaa  are  llluetrated  in  Fig.  8.  The  third  paaa  produced  aimilar  reaulta;  the  firat  paae  waa  too  poor 
to  be  of  uae.  Returne  from  the  Culf  of  St.  Lawrence  were  not  analyaed, 

4.2  April  8,  1982 

Four  paaaea  were  made  over  the  coverage  area  during  the  aix  houra  of  obaervation  on  April  8, 
1982.  Uae  of  the  re-deaigned  array  configuration  provided  coverage  of  the  entire  area  without 
amblguoua  reeponaea. 

The  wave  height  reaulta  were  aparae.  All  of  the  acceptable  reaulta,  from  all  four  obaervation 
paaaea,  are  ehown  in  Fig.  9,  overlaid  on  a hindcaet  wave-height  contour  map,  The  radar  returna  were 
recorded  ovar  the  period  14:002  to  20:002;  the  contour  map  waa  compiled  effective  12:002.  The 
Atlantic  atom,  with  11  metre  wavea  at  ita  centre,  waa  approaching  from  tha  aouth.  The  radar 
maaeuramenla  agreed  quite  well  wtth  the  hindcaet  map.  There  waa  only  one  conapicuoualy  high  value. 

Tha  wave  heighta  e-own  within  the  lend  area  of  northern  Newfoundland,  which  are  thought  to  be  the 
reault  of  aidelobe  reeponaea,  were  in  agreement  with  the  low  wave  heighta  near  the  tip  of  the  leland. 

Tha  corraapondlng  wave-direction  map,  Fig,  10,  ahowa  good  agreement  with  expected  wind 
directiona.  The  preaaure  iaobara  were  plotted  for  an  effective  time  of  18:002,  At  waa  the  caae  for 
wave  heighta,  wave  directiona  apparently  raeaeured  over  land  arete  agreed  with  wind  directiona  expected 
over  nearby  ocean  areat:  The  wave  direction  mtaauramanta  produced  more  than  one  reault  for  moat 

meaeurement  cella,  but  tinea  thaaa  all  ahowed  general  agreamant,  only  one  waa  plotted  for  each  cell. 

Aaaaaament  of  radar  errora,  by  compariaon  of  radar  raaulta  againat  meteorological  data  and 
chip  reporta,  ia  made  difficult  by  uncertainty  in  tha  turface  data  and  by  difference  in  the  time  of 
obaervation.  For  example,  tome  of  the  wave  height  reaulta  plotted  in  Fig.  9 were  recorded  more  than  6 
houra  attar  tha  effective  time  of  the  hindcaet  wave  height  contoura,  When  that  data  waa  extracted  and 
compared  againat  the  more  timely  hindcaat  map  compiled  for  24:002,  it  waa  found  that  the  agreement  waa 
poor  at  the  aouthern  limit  of  the  area.  Tho  aurface  data  indicated  that  tha  approaching  atorm  had 
entered  the  area,  while  the  radar  data  indicated  that  it  had  not,  Some  qualification  ia  neceaaary 
whenever  conpariaona  are  made. 


5,0  RESULTS,  C-CORE 


The  map  ahown  in  Fig.  11  ia  a compoalte  of  accepted  wave  directiona  from  aalected  areae  of 
opan  ocean  obaerved  on  March  26,  1982.  Extenaive  ice  cover  along  the  coaatline  reduced  the  number  of 
obaervation  cella  to  81,  from  which  47  yielded  wave  directiona.  The  inherent  right/laft  ambiguity  of 
the  radar  raaulta  waa  reaolvad  by  compariaon  to  meteorological  data.  Only  one  wave  direction  waa 
plotted  for  each  cell,  although  up  to  three  reaulta  were  obtained  from  analyaia  of  the  three  radar 
paaaea  over  the  area.  In  all  but  three  cella  the  reaulta  ware  coneiatant. 

The  difference  in  wave  direction  diaplayed  by  the  two  aoutherniroat  obaervation  areaa  waa 
coneiatant  with  an  aroa  of  high  atmoapheric  preaaure  located  juat  off  the  coutheaatern  corner  of  the 
map, 

Thia  map  waa  produced  from  what  ia  thought  to  be  one  of  the  beat  obaervinj  daya.  Even  ao, 
more  than  three-quartere  of  the  apectra  produced  were  coneidered  to  be  unacceptable.  From  tne  327 
apectra  obtained  from  the  82  target  cella,  only  79  were  accepted.  Moat  of  the  rejected  apectra 
diaplayed  either  aignificant  lonoapheric  contamination  or  inadequate  aignal  atrength.  In  acme  caaea 
it  waa  obvioua  that  the  apectra  bad  not  been  correctly  aligned  and  conaequently  the  Bragg  peak  had  not 
been  correctly  identified. 

Some  of  the  analyaia  problema  are  llluetrated  by  the  unrelated  apectra  ahown  in  Fig.  12. 

Although  apectra  A and  B are  quite  aimilar,  both  dlaplaying  readily  identifiable  Bragg  peaka,  apectrum 
B waa  miaaligned  by  failure  of  the  centering  algorithm,  which  incorrectly  identified  the  right-hand 
limit  at  the  point  Indicated  by  the  arrow.  Spectra  C and  0 both  diaplayed  multipath  contamination  and 
eo  were  rejected  after  viaual  inapection. 

It  waa  realiaed  that  incluaion  of  all  initial  apectra  into  the  averaging  proceaa  would  reault 
in  a high  rejection  rate  of  the  final  averagea,  but  it  waa  felt  that  the  uae  of  a ahort  obaervation  >1 

time  and  the  incluaion  of  only  a araall  number  of  range  cella  would  uaually  provide  an  acceptable  final  j 

average  apectrum.  However,  it  appeara  that  the  rajection  rate  may  be  etill  too  high  to  build  adequate  1 

wind  direction  mapa  from  the  available  data.  1 
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FIG.  9.  Wave  Heights  Measured  on  April  8.  19 B2 

At  higher  wave  heights,  the  radar  results  agree  more  closely  with  the  meteorological  data. 
The  radar  data  wao  recorded  between  IAiOOZ  and  20 : 00Z;  the  hindcaat  wave  height  contours 
were  plotted  effective  12:00Z. 
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FIG.  10,  wave  Directions  Measured  on  April  6,  1982 


The  radar  data  were  recorded  between  14iG0Z  and  20:002}  the  hindcaat  meteorological 
isobars  were  plotted  effective  18:002, 
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FIG.  111  Cgagogltg  Wav*  D1 ruction  Map,  March  26.  1982 

Th*  arrow*  lndlcata  tha  wav*  direction*  for  result*  retained  after  viaual  inspection  of 
the  apaotra.  Each  arrow  indicate*  tha  on*  raault  obtained  from  th*  averaging  of  th* 
apectra  from  8 adjacent  range  cal la . Th*  alee  of  th*  final  target  cell  la  indicated  by 
the  dotted  area*  ahown  on  tha  vaatam  and  of  tha  candle  area*. 
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FIG.  12.  Sony  1*  a of  Spaotra  from  Automatic 

Frocaaalng  of  On*  of  th*  Bettor-Quality 
Data  Sate,  March  26,  1982 

Th*  vertical  line*  lndlcata  tha  poaltlona  of  tha 
Bragg  peak*  and  tha  cantrold  of  th*  apactrum  aa 
identified  by  tha  algorithm*.  Th*  horlcontal  line 
lndlcata*  tha  average  value  (power)  of  th*  apactrum, 
A proposed  quality  index  la  reported  aa  "PKD", 
higher  value*  indicating  batter  quality. 

From  these  examples,  apactrum  A was  tha  only  one 
accepted.  Although  apeerrum  B was  eladlar  In 
quality,  tha  search  strategy  failed  because  of 
th*  low  point  indicated  by  tha  arrow. 
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6.0  CONCLUSIONS 


At  a «kyvav«  itt-itici  radar  workehop  held  In  Rockville,  Maryland,  in  May,  1981  [IS],  required 
capabilitiaa  (or  operational  akywave  radara  vere  propoeed.  Among  thaae  were: 


(i)  accuracy  of  meaauranent : 


(ii) 


wave  height: 
wave  direction: 
poaition  fixeai 
frequency  of  reporta : 


M).5m  or  ^10t 
+ 20  degree* 

+ 25  km 
twice  dally. 


From  the  limited  amount  of  data  ao  far  analyaed  in  thie  experiment,  there  la  eome  baaia  for 
optimiam  that  the  accuracy  requirement#  for  wave  height  and  wave  direction  may  be  met,  even  at  tha 
geomagnetic  latitudaa  of  the  Canadian  environment,  Meaeurement  of  wave  direction  doee  not  appear  to 
be  a problem}  the  manual  technique  derived  vat  id  reaulta  even  from  unconditionel  average#  of  the 
apactra,  and  an  automatic  analyaia  technique  appear#  to  be  feaelble.  Accurate  meaeurement  of  low  wave 
height#  (2  or  3 metre#)  haa  been  found  to  ba  difficult  becauae  of  tha  aenaltivity  of  tha  procaae  to 
multipath  contamination. 


The  problem  of  accurate  poaition  flxaa,  under  atudy  but  not  reported  hare,  may  preeent  eome 
difficult!##.  Poaition  fixaa  depend  upon  accuracy  in  the  determination  of  both  the  beam  direction  and 
the  virtual  height  of  the  ionoaphara.  Direction-finding  experience  with  lonoapheric  tilte,  and  with 
remote  eatlmataa  of  virtual  height,  auggaata  that  the  +23  km  apecification  may  not  ba  attainable  for 
radar  aignala  propagated  via  tha  F region.  ” 

The  moat  formidable  problem,  however,  ia  likely  to  be  the  requirement  for  twice-daily 
coverage.  Aa  indicated  earlier,  tha  reaulta  given  hare  ware  aelectad  becauae  they  were  aucceaaful. 
Although  tha  aaopa  of  thia  experiment  waa  not  aufficiant  to  yield  an  aatimata  of  reliability  (even 
with  more  complete  analyaia),  tha  problem#  of  finding  good  operating  daya  during  the  execution  of  the 
experiment  indicate  that  there  would  ba  aignificant  gapa  in  coverage  by  an  operational  radar. 
Night-time  operationa  have  not  bean  examined  at  all.  However,  the  U.S.  worker#,  at  the  Wide  Aperture 
Receiving  Facility  in  California  (16),  have  demonatrated  a capability  for  more-or-leea  routine 
operation,  at  1 -eat  at  their  geomagnetic  latitude#.  Their  technique  ia  dependent  upon  a capability 
for  on-line  real-time  analyaia,  which  permit#  tha  experimenter  to  peralat  until  he  aucceeda,  and  upon 
freedom  of  choice  of  operating  frequency. 

It  ia  hoped  that  continued  effort  to  analyae  the  remainder  of  the  data  recorded  in  thia 
experiment  will  ahed  more  light  on  all  of  theae  queatlnna. 
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